Bioaccumulation modelling and sensitivity analysis for discovering key players in contaminated food webs: the case study of PCBs in the Adriatic Sea

Modelling bioaccumulation processes at the food web level is the main step to analyse the effects of pollutants at the global ecosystem level. A crucial question is understanding which species play a key role in the trophic transfer of contaminants to disclose the contribution of feeding linkages and the importance of trophic dependencies in bioaccumulation dynamics. In this work we present a computational framework to model the bioaccumulation of organic chemicals in aquatic food webs, and to discover key species in polluted ecosystems. As a result, we reconstruct the first PCBs bioaccumulation model of the Adriatic food web, estimated after an extensive review of published concentration data. We define a novel index aimed to identify the key species in contaminated networks, Sensitivity Centrality, and based on sensitivity analysis. The index is computed from a dynamic ODE model parametrised from the estimated PCBs bioaccumulation model and compared with a set of established trophic indices of centrality. Results evidence the occurrence of PCBs biomagnification in the Adriatic food web, and highlight the dependence of bioaccumulation on trophic dynamics and external factors like fishing activity. We demonstrate the effectiveness of the introduced Sensitivity Centrality in identifying the set of species with the highest impact on the total contaminant flows and on the efficiency of contaminant transport within the food web

Estimating dormant and active hematopoietic stem cell kinetics through extensive modeling of bromodeoxyuridine label-retaining cell dynamics.

Bone marrow hematopoietic stem cells (HSCs) are responsible for both lifelong daily maintenance of all blood cells and for repair after cell loss. Until recently the cellular mechanisms by which HSCs accomplish these two very different tasks remained an open question. Biological evidence has now been found for the existence of two related mouse HSC populations. First, a dormant HSC (d-HSC) population which harbors the highest self-renewal potential of all blood cells but is only induced into active self-renewal in response to hematopoietic stress. And second, an active HSC (a-HSC) subset that by and large produces the progenitors and mature cells required for maintenance of day-to-day hematopoiesis. Here we present computational analyses further supporting the d-HSC concept through extensive modeling of experimental DNA label-retaining cell (LRC) data. Our conclusion that the presence of a slowly dividing subpopulation of HSCs is the most likely explanation (amongst the various possible causes including stochastic cellular variation) of the observed long term Bromodeoxyuridine (BrdU) retention, is confirmed by the deterministic and stochastic models presented here. Moreover, modeling both HSC BrdU uptake and dilution in three stages and careful treatment of the BrdU detection sensitivity permitted improved estimates of HSC turnover rates. This analysis predicts that d-HSCs cycle about once every 149-193 days and a-HSCs about once every 28-36 days. We further predict that, using LRC assays, a 75%-92.5% purification of d-HSCs can be achieved after 59-130 days of chase. Interestingly, the d-HSC proportion is now estimated to be around 30-45% of total HSCs - more than twice that of our previous estimate

Trends in modeling Biomedical Complex Systems

In this paper we provide an introduction to the techniques for multi-scale complex biological systems, from the single bio-molecule to the cell, combining theoretical modeling, experiments, informatics tools and technologies suitable for biological and biomedical research, which are becoming increasingly multidisciplinary, multidimensional and information-driven. The most important concepts on mathematical modeling methodologies and statistical inference, bioinformatics and standards tools to investigate complex biomedical systems are discussed and the prominent literature useful to both the practitioner and the theoretician are presented

An integrated modelling approach for R5-X4 mutation and HAART therapy assessment

We have modelled the within-patient evolutionary process during HIV infection using different methodologies. New viral strains arise during the course of HIV infection. These multiple strains of the virus are able to use different coreceptors, in particular the CCR5 and the CXCR4 (R5 and X4 phenotypes, respectively)influence the progression of the disease to the AIDS phase. We present a model of HIV early infection and CTLs response which describes the dynamics of R5 quasispecies, specifying the R5 to X4 switch and effects of immune response. We illustrate dynamics of HIV multiple strains in the presence of multidrug HAART therapy. The HAART combined with X4 strain blocker drugs might help to reduce infectivity and lead to slower progression of disease. On the methodology side, our model represents a paradigm of integrating formal methods and mathematical models as a general framework to study HIV multiple strains during disease progression, and will inch towards providing help in selecting among vaccines and drug therapies. The results presented here are one of the rare cases of methodological cross comparison (stochastic and deterministic) and a novel implementation of model checking in therapy validation

Modelling circulating tumour cells for personalised survival prediction in metastatic breast cancer.

Ductal carcinoma is one of the most common cancers among women, and the main cause of death is the formation of metastases. The development of metastases is caused by cancer cells that migrate from the primary tumour site (the mammary duct) through the blood vessels and extravasating they initiate metastasis. Here, we propose a multi-compartment model which mimics the dynamics of tumoural cells in the mammary duct, in the circulatory system and in the bone. Through a branching process model, we describe the relation between the survival times and the four markers mainly involved in metastatic breast cancer (EPCAM, CD47, CD44 and MET). In particular, the model takes into account the gene expression profile of circulating tumour cells to predict personalised survival probability. We also include the administration of drugs as bisphosphonates, which reduce the formation of circulating tumour cells and their survival in the blood vessels, in order to analyse the dynamic changes induced by the therapy. We analyse the effects of circulating tumour cells on the progression of the disease providing a quantitative measure of the cell driver mutations needed for invading the bone tissue. Our model allows to design intervention scenarios that alter the patient-specific survival probability by modifying the populations of circulating tumour cells and it could be extended to other cancer metastasis dynamics